Rashba spin-orbit-interaction-based quantum pump in graphene 
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We present a proposal for an adiabatic quantum pump based on a graphene monolayer patterned 
by electrostatic gates and operated in the low-energy Dirac regime. The setup under investigation 
works in the presence of inhomogeneous spin-orbit interactions of intrinsic- and Rashba-type and 
allows to generate spin polarized coherent current. A local spin polarized current is induced by the 
pumping mechanism assisted by the spin-double refraction phenomenon. 
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Graphene-based devices have recently been the ob- 
ject of intense experimental and theoretical studies. 
Due to its highly integrable ultra-flat geometry, 1 high 
electron mobility, tunable carrier concentration, 2 high 
chemical homogeneity, and excellent intrinsic transport 
properties, 3 graphene monolayers are becoming the key 
technological ingredient in molecular electronics. On the 
other hand, quasi-particles in the honeycomb lattice of 
graphene — close to the charge neutrality point — can be 
described as massless Dirac- Weyl fermions. All these ex- 
citing properties make graphene a promising material for 
realizing devices having exotic functionalities. An impor- 
tant class of these devices are the pure spin-current gen- 
erators (PSGs), being them conceptual relevant devices 
in semiconductor spintronics. 4 Among different spintron- 
ics mechanisms that are useful for the design of a PSG, 
an interesting option is the implementation of the adi- 
abatic quantum pumping (AQP) technique. 5 However, 
implementation of PSG schemes is strongly material- 
dependent and therefore requires an accurate control 
of interface properties. Motivated by this requirement, 
the application of the AQP to graphene-based devices 
has become a relevant research activity in spintronics- 
oriented molecular electronics. 6 In addition, the impor- 
tant role of the spin-orbit interaction (SOI) and its tun- 
ability via all-electrical means in graphene-based systems 
has recently been reported. 7 Despite the relevance of the 
SOI for spintronics, the latter has not yet been intensely 
studied in combination with the AQP technique in order 
to provide a mechanism to generate a pure spin current. 

Here, we present a proposal of AQP for graphene in 
the presence of gate tunable SOI — we show this setup 
may lead to spin-polarized ballistic transport. 

Our proposal is based on the exploitation of the spin- 
double refraction phenomenon 8 (SDRP), which occurs 
at the interface between the SOI-active and the normal- 
graphene regions, in the AQP framework. Other proposal 
of quantum pumping in graphene recently appeared 9-11 
but none of these considered the relevant effect of control- 
lable SOIs. Symmetries in graphene allow for two types 
of SOIs: 12 the intrinsic and the Rashba one. The former, 
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FIG. 1: Graphene-based pumping device. The monolayer 
graphene is sandwiched between two metallic contacts. An 
electrostatic top gate (region T) of length Lt and a back gate 
(region SO) of length Lso are placed between the two con- 
tacts, which leave two 'normal' graphene regions Ni and Na- 
Finite spin-orbit interactions are present in region SO, only. 
The top-part of the figure schematically shows the relevant 
kinematic angles fa for the scattering process. In the SO 
region, spin double refraction is observed, i.e., the two spin 
species move in different directions (f>f. 



originating from carbon intra-atomic SOI, has been esti- 
mated to be rather weak in clean flat graphene. 13 How- 
ever, a recent theoretical proposal suggests how atoms- 
mediated hopping can enhance it. 14 The Rashba SOI 
originates from interactions with the substrate, Stark 
and/or curvature effects. 13 It has been experimentally 
enhanced in graphene samples on Ni with intercalated 
Au atoms. 7 

We consider a graphene-based system in the low- 
energy regime assuming that the two Dirac cones can be 
treated independently, meaning that they are not mixed 
by interaction effects or boundary conditions. 15 Within 
this approximation the single-valley Hamiltonian for the 
system reads: 

X(x) 

H =v F (s (g>a) -p-\ — (s y ®<t x -8 x ® ay) 

+ A{x){s z ®o- z )+V T (x) + Vso{x). (1) 
Here, the cr x / y / z and s x / y / z are Pauli matrices describing 
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the pseudo-spin and fermionic spin degree-of-freedom, re- 
spectively. The matrix ctq — sq — I2 is the identity ma- 
trix of dimension 2. In the Hamiltonian (1) the first term 
describes particles with Dirac-like energy dispersion, the 
second one is the Rashba SOI, while the third term takes 
into account the intrinsic SOI of graphene. Let us note 
that at the lowest order in the long-wavelength approx- 
imation, the velocity operator v = v-p(so Cg> a) does not 
involve spatial derivatives. 8 ' 16 The fourth term in (1) de- 
scribes an electrostatic gate of length Lt and strength Vr 
placed before the region of finite SOI. The latter term de- 
scribes a region of length L$o characterized by the pres- 
ence of a back-gate V30 having the two-fold action of (i) 
varying the electrochemical potential by a quantity Vso, 
and (ii) tuning the intensity of the Rashba SOI according 
to X(x) « Xq(x) + SgVso(x). Here Sg is a proportional- 
ity constant between the Rashba SOI and the applied 
electrostatic potential Vso, while Ao is the bare Rashba 
SOI. 13 

In order to study AQP in this set-up, we allow the 
electrostatic potentials in (1) to be time-dependent quan- 
tities. In particular, their time-dependence is chosen as 
Vf){x,t) = V${x) + 6Vpsm(ut + ip p ) with G {T,SO}, 
where ip = ipso — 7^ is a phase difference between the 
two pumping parameters. An important consequence of 
the gate modulation is an induced time-dependence also 
for the Rashba SOI (X(x) X(x,t)). 

Assuming current conservation laws, 17 the charge cur- 
rent density Jq and the spin current density J§ can be 
expressed as Jq = qv-p^ (s <8> (r)tp and Js = ^v-p^ (s z ® 
respectively. In order to find an expression for the 
pumped charge and spin currents, we evaluate the scat- 
tering matrix S of the whole system. Thus, we require 
continuity of the particle wave-function at all potential 
discontinuities of Eq. (1) along the x axis. We consider 
hereafter a purely electronic transport in the leads, while 
the SOI parameters (both A and A) are taken finite only 
in the scattering region labeled 'SO' in Fig. 1. Charge 
and spin currents generated in the system can be evalu- 
ated within a modified scattering field approach. 20,21 The 
currents in lead a = 1, 2 are obtained by the average over 
incoming angles (pi ni 
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7 = C/S, and the </>; n -dependent current contributions 
read 



J" {<Pin) = (2A/1 sin(<p) cos 
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where the scattering matrix is also a function of (pin- 
Here, a factor of two accounts for the valley degener- 
acy and Af± = (EpW) / '(hvp) is the number of transverse 
channels. The constant A^ takes the value A^ = q for the 
charge current and Af = sh/2 for the spin current. Note 



that the pumped currents are maximal for a phase differ- 
ence ip = it/2 in the time variation of the two pumping 
parameters Vt and Vgo • In the following we will express 
the currents in units of ^ A/j_ ^iP , with a pre- factor q 
for the charge current and H/2 for the spin currents. The 
order of magnitude of the pumped charge current is of the 
order 1-20 nA assuming w/(2ir) = 5 GHz, Af± = 100 and 
the weak pumping limit SVp/Vp ~ 0.1. For the Rashba 
SOI strength we take the maximal value of A = 13 meV. 7 
We start by analyzing the refraction processes at 
each of the three interfaces Ni-T, T-SO, and SO-N 2 
(c.f. Fig. 1). Translational invariance along the y di- 
rection implies conservation of the momentum along the 
interfaces, i.e., ^fe erface ) = k y (x+ tcrlacc ). By using 
simple kinematic considerations 8 we can write a general 
relation expressing the refraction angle at the interfaces. 
This reads 



_l = arcsm 
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where the fc$ and <pi (with i = 1,...,4 and cp\ = <f>- m ) 
are the momentum modulus and propagation directions 
in the four regions, respectively (c.f. Fig. 1). At each 
interface there exists an incoming critical angle (pi over 
which the incoming wave is totally refracted along the 
interface, this is given by 
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(5) 



When entering a region of finite Rashba SOI, an in- 
coming wave splits into two outgoing waves with differ- 
ent propagation directions which are associated to the 
two different spin eigenstates within the Rashba SOI re- 
gion. This refraction along the different spin channels is 
called SDRP 8 and is relevant for the spin polarized trans- 
port. When one of the modes reaches the critical angle, 
it propagates completely parallel to the interface, while 
becomes evanescent in the perpendicular direction. 8 

We present results for two cases: the first one with 
negligible intrinsic SOI and a second one with a finite 
intrinsic SOI, still smaller than the Rashba SOI. In gen- 
eral the spin [charge] pumping current is an odd [even] 
function of the injection angle </>i n , as can be seen from 
the symmetry properties of the SOI region. 8 ' 18 

In Fig. 2 we report the angular dependence of the 
charge and spin pumping currents for different lengths 
of the SOI region. For small values of the injection an- 
gle, both currents are rather small, then for > 0.5 
they both acquire finite values and oscillate as a function 
of the injection angle. For injection angles greater than 
the first critical angle of the Ni-T interface they become 
negative with a peak at the smaller of the two critical 
angles </>2 ,=t of the T-SOI interface. For even larger val- 
ues of 4>i n , the two currents both approach zero. A finite 
value of the intrinsic SOI increases the overall value of 
the currents. 
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FIG. 2: Charge Jc (left panels) and spin J7s (right panels) 
pumping currents as a function of the injection angle 0i n = (pi 
for different lengths of the SOI region Lso: (i) 50 nm (solid 
line), (ii) 75 nm (dashed line), (hi) 100 nm (dotted-dashed 
line). The others parameters are: = 1.3 meV, Lt = 
50 nm, E in = 71.5 meV, A = 11.7 meV, A ~ meV (upper 
panels) and A ~ 2.6 meV (lower panels), Vg°o = 13 meV, and 
6g = 0.1. 



average charge pumping current is zero while the spin 
current is finite. For the case of a 100 nm SOI region 
[panels (c) and (f)] we can identify two of these working 
points. Furthermore, we observe that a finite value of the 
intrinsic SOI is producing only a change in the magnitude 
of the currents. 

In conclusion, we have studied a graphene-based spin- 
tronic pump exploiting the SDRP. This is due to the 
inhomogeneous Rashba SOI that can be obtained via an 
appropriate spatial modulation of the Au atoms between 
the graphene layer and the Ni substrate. 7 Special signa- 
tures of the critical refraction angles have been identified 
in the angle-resolved properties of the pump. Charge 
and spin current curves as a function of the gate volt- 
ages support the possibility to obtain local spin-polarized 
transport in the absence of a charge current. An imaging 
technique similar to the one described by Topinka et al. 22 
can be used to detect the coherent spin-polarized particle 
flow with angular resolution. 

We acknowledge N. Andrei, A. Di Bartolomeo and L. 
Lenz for useful discussions. The work of DB is supported 
by the Excellence Initiative of the German Federal and 
State Governments. 



The regions where the charge current changes sign are 
particularly interesting. There, by performing an inte- 
gration of both currents over a small range of injection 
angles, it is possible to find a considerable value for the 
pumped spin current while the pumped charge current 
is negligible. The evaluation of the average currents by 
integration over a small range of angles can be experi- 
mentally implemented by using a quantum point contact 
(QPC) far away from the scattering region. By chang- 
ing the width of the QPC, which changes the angular 
structure of the electron flow, it is possible to change 
the angular contribution to the transmission through the 
contact itself. 22 The data shown in Fig. 2 allows to iden- 
tify distinct regions with this property to be considered 
as working points of a PSG (c.f. caption of Fig. 3 for their 
actual value). 

The average pumping currents are shown in Fig. 3 as 
a function of the electrostatic voltage Vso for different 
lengths of the SOI region with and without intrinsic SOI. 
We recall the twofold action of Vso which produces a 
chemical potential shift accompanied by a variation of 
the Rashba SOI strength. The main result is an increase 
of the overall currents for an increase of the SOI region 
length. In addition, we observe the existence of optimal 
working points in terms of Vso for the PSG in which the 
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FIG. 3: Angle averaged charge (J7c) (blue-solid line) and spin 
(c7s) (red-dashed line) pumping currents as a function of the 
electrostatic gate voltage Vso f° r different lengths I/so of the 
SOI region: (a) and (d) 50 nm, (b) and (e) 75 nm, and (c) and 
(f ) 100 nm. The integration region is </>j n 6 [0.4, 0.57] for (a) 
and (d), <j) in G [0.61,0.69] for (b) and (e) and <j> in G [0.9, 1.05] 
for (c) and (f). For the upper panels (a)-(c) A = 0, for the 
lower ones (d)-(f) A = 2.6 meV. The other parameters are 
common to all the panels: Vx = 1-2 meV, Lt = 50 nm, 
E in = 71.5 meV, A = 11.7 meV, and 5g = 0.1. 
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